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[57] ABSTRACT 

A method for maintaining the attitude of a three -axis 
controlled satellite by use of magnetic torquers includes 
using magnetometers for measuring the direction of the 
ambient geomagnetic field. The direction of the net 
reaction wheel momentum is also determined. The 
angle between the direction of the geomagnetic field 
and the net reaction wheel momentum is determined. 
The angle is compared with a threshold value. Mag- 
netic torquer power consumption is reduced by operat- 
ing the magnetic torquers only when the angle exceeds 
the threshold value. 

5 Claims, 7 Drawing Sheets 
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SPACECRAFT MOMENTUM UNLOADING USING 
CONTROLLED MAGNETIC TORQUES 

The invention described herein was made in the per- 
formance of work under NASA Contract No. NAS5- 
32000 and is subject to the provisions of Section 305 of 
the National Aeronautics and Space Act of 1958, as 
amended (72 STAT. 435; 42 U.S.C. 2457). 

BACKGROUND OF THE INVENTION 

This invention relates to attitude control of space- 
craft, and more particularly to unloading of reaction or 
momentum wheel momentum by the use of controlled 
magnetic torquers. 

Orbiting spacecraft are used for a large variety of 
sensing and communication purposes for photographic 
purposes, it is desirable for the spacecraft to be rela- 
tively near the Earth, so that the cameras or sensors are 
close to the subjects, and so that the spacecraft may 
cover large amounts of surface by virtue of a rapid 
orbital rate. For communication purposes, a geosyn- 
chronous equatorial orbit is often desirable. Whatever 
the orbit, the satellite must be stabilized in space if the 
sensors or antennas are to be pointed in appropriate 
direction. 

Spacecraft attitude stabilization may be accomplished 
by spinning the spacecraft and by mounting the sensors 
or antennas on a despun platform. Alternatively, the 
spacecraft may be stabilized in three axes. Three-axis 
stabilization may be accomplished by a control system 
using fuel-burning thrusters, but the use of such thrust- 
ers requires the expenditure of fuel, which tends to limit 
the service life of the spacecraft. Another method of 
three-axis stabilization uses reaction wheels or momen- 
tum wheels mounted within the spacecraft. The reac- 
tion wheels are torqued under the control of an attitude 
control system to cause the spacecraft body to assume 
the desired attitude. The energy for torquing the wheels 
is electrical, and is derived from solar panels. Attitude 
control by the use of wheels therefore generally pro- 
vides longer operating lifetime than simple attitude 
control by the use of thrusters. 

In order to maintain a face of the satellite directed 
toward the Earth, the satellite must rotate about its 
pitch axis at the orbit rate. When using reaction wheel 
control, the wheel or wheels are torqued to maintain an 
earth-oriented attitude. Environmental disturbance 
torques that act on the spacecraft tend to result in a 
gradual increase in the wheel momenta. In order to 
prevent the wheels from reaching their mechanical 
momentum limits, and to thereby maintain the ability to 
provide reaction wheel based attitude control, the 
wheel momentum must be periodically reduced. Wheel 
momentum may be reduced by the use of fuel-consum- 
ing thrusters which torque the spacecraft body in a 
manner which opposes the wheel momentum. How- 
ever, the use of such thrusters is subject to the same 
objection of limited lifetime described above. 

FIG. la illustrates a spacecraft designated generally 
as 10, and including a body 12 and a solar panel 14. 
Spacecraft 10 is following an orbit illustrated by 
dashed-line path 16 about the Earth 18. The Earth’s 
rotational axis is illustrated as 20, and the Equator is 
illustrated by dash-line 22. Lines of magnetic force asso- 
ciated with the Earth’s magnetic field are illustrated as 
24. As illustrated in FIG. 1, solar panel 14 is directed 
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toward the sun 26, and is controlled to continue facing 
the sun at all orbital positions of satellite 10 . 

FIG. lb illustrates details of spacecraft 10. In FIG. lb 
a first reaction wheel 30 is oriented with its axis along 
5 the pitch or y axis, a second reaction wheel 32 is ori- 
ented with its axis along the x or roll axis, and a third 
reaction wheel 34 is oriented with its axis along the yaw 
or z axis. Also in FIG. lb, magnetic torquer windings 
illustrated as 36, 38 and 40 are wound about the space- 
10 craft body. Further in FIG. lb, a set of three-axis mag- 
netometers is illustrated as a block 41, and is connected 
to a torquer control arrangement illustrated as a block 
42 

FIGS, la, lb and 2c illustrate as time functions the 
15 roll, pitch and yaw disturbance torques, respectively, 
measured in inch-lb. In FIG. 2, the disturbance torques 
include a periodic component and an offset component. 
The disturbance torques arise from atmospheric drag, 
gravity gradient effects, and solar pressure. The cyclical 
20 perturbations recur at roughly 6,000 second intervals, 
corresponding to the orbital period. The constant offset 
of the disturbance torques is attributable to the asym- 
metric orientation of solar panel 14 of FIG. la. In FIG. 
2, zero seconds and each 6,000 second increment there- 
25 after represent the times at which the spacecraft crosses 
the Earth’s Equator, known as the ascending node. 

FIGS. 3 a, 3 b and 3c represent as time functions the 
roll, pitch and yaw wheel momentum, respectively, in 
the absence of unloading. The changes in wheel mo- 
30 mentum with time are attributable to the disturbance 
torques illustrated in FIG. 2. Some components of the 
disturbance torques of FIG. 2 contribute to the cyclical 
variations of wheel momentum in FIG. 3, while others 
result in secular or continually increasing or decreasing 
35 components. For example, in FIG. 3c, the peak values 
of the cyclical component of wheel momentum in- 
creases with increasing time. At some point, therefore, 
such a continuing increase in amplitude will result in 
reaching the wheel momentum limit. Similarly, the 
40 pitch wheel momentum of FIG. 3b increases monotoni- 
cally with increasing time. In order to maintain continu- 
ous attitude control, the secular component of the 
wheel momentum must be unloaded. 

FIGS. 4 a, 4 b and 4c represent the roll, pitch and yaw 
45 torquer magnetic moment measured in ampere-tums- 
meter 2 (ATM 2 ) for a spacecraft using a prior art torquer 
control system which pulse width or continuously mod- 
ulates the torquer current. The prior art control system 
measures the Earth’s geomagnetic field B and the net 
50 momentum of the reaction wheel, and forms a vector 
cross product HXB. The normalized HXB vector is 
multiplied by the maximum magnetic moment which 
can be achieved along each of the roll, pitch and yaw 
axes to determine the torquer dipole command. In 
55 FIGS. 4 a, 4 b and 4c, large excursions or jitter in the roll 
and pitch moments, and lesser excursions in the yaw 
moments, occur at about 7500 seconds, 14,000 seconds, 
and 20,000 seconds, corresponding to passage of the 
spacecraft across the Earth’s magnetic poles. The jitter 
60 illustrated in FIG. 4 arises because the magnitude of the 
Earth’s geomagnetic field near the poles is large, but the 
cyclical component of the wheel momentum results in 
relatively low wheel momentum. Operation of the mag- 
netic torquers in regions near the poles, therefore, pro- 
65 duces a high magnitude of wheel momentum change at 
times when the wheel momentum is low, which tends to 
drive the torquer command from one extreme to the 
other. This in turn results in the illustrated jitter, which 
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disturbs the spacecraft attitude. Such disturbances are 
undesirable. 

SUMMARY OF THE INVENTION 

A method for maintaining the attitude of the satellite 
includes the step of measuring the direction of the ambi- 
ent magnetic field in the vicinity of the satellite. The 
magnitude and direction of the net wheel momentum of 
the satellite is also determined, and the angle therebe- 
tween is determined. The angle is compared with a 
predetermined value, and magnetic torquers are ener- 
gized only when the angle exceeds a predetermined 
value. In a particular embodiment of the invention, the 
magnitude of the ambient magnetic field is also mea- 
sured. 

DESCRIPTION OF THE DRAWING 

FIG. la is a perspective or isometric view of an asym- 
metric satellite in a polar orbit around the Earth, with 
its solar panel pointed toward the Sun, and FIG. 16 
illustrates the orientation of reaction or momentum 
wheels and magnetic torquing coils in the spacecraft of 
FIG. lb, FIGS, la and lb are referred to jointly as FIG. 

i; 

FIGS. 2a, lb and 2c, referred to jointly as FIG. 2, are 
plots versus time of environmental disturbance torques 
about roll, pitch and yaw axes, respectively; 

FIGS. 3a, 3 b and 3c, referred to jointly as FIG. 3, are 
plots versus time of wheel momentum about roll, pitch 
and yaw axes, respectively, without momentum unload- 
ing control; 

FIGS. 4a, 4 b and 4c, referred to jointly as FIG. 4, are 
plots versus time of magnetic moments produced about 
roll, pitch and yaw axes, respectively, by spacecraft 
torquers using a prior-art control system; 

FIG. 5 is a simplified block diagram of a control 
system according to the invention; 

FIGS. 6, 7, 8 and 9 a are simplified flow charts illus- 
trating details of the control system of FIG. 5, and FIG. 
9b illustrates various angles associated with the flow 
chart of FIG. 9a; 

FIGS. 10a, 10 b and 10c, referred to jointly as FIG. 

10, are plots versus time of magnetic moments produced 
about roll, pitch and yaw axes, respectively, by space- 
craft torquers using a control system according to the 
invention; 

FIGS. 11a, 116 and 11c, referred to jointly as FIG. 

11 , are plots versus time of wheel .momentum about roll, 
pitch and yaw axes, respectively, using momentum 
control according to the invention. 

DESCRIPTION OF THE INVENTION 

FIG. la illustrates a spacecraft designated generally 
as 10, and including a body 12 and a solar panel 14. 
Spacecraft 10 is following an orbit illustrated by 
dashed-line path 16 about the Earth 18. The Earth’s 
rotational axis is illustrated as 20, and the Equator is 
illustrated by dash-line 22. Lines of magnetic force asso- 
ciated with the Earth’s magnetic field are illustrated as 
24. As illustrated in FIG. 1, solar panel 14 is directed 
toward the sun 26, and is controlled to continue facing 
the sun at all orbital positions of satellite 10 . 

FIG. 16 illustrates details of spacecraft 10. In FIG. 16 
a first reaction wheel 30 is oriented with its axis along 
the pitch or y axis, a second reaction wheel 32 is ori- 
ented with its axis along the x or roll axis, and a third 
reaction wheel 34 is oriented with its axis along the yaw 
or z axis. Also in FIG. 16, magnetic torquer windings 
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illustrated as 36, 38 and 40 are wound about the space- 
craft body. Further in FIG. 16, a set of three-axis mag- 
netometers is illustrated as a block 41, and is connected 
to a torquer control arrangement illustrated as a block 
5 42 

FIGS. 2a, 26 and 2c illustrate as time functions the 
roll, pitch and yaw disturbance torques, respectively, 
measured in inch-lb. In FIG. 2, the disturbance torques 
include a periodic component and an offset component. 
10 The disturbance torques arise from atmospheric drag, 
gravity gradient effects, and solar pressure. The cyclical 
perturbations recur at roughly 6,000 second intervals, 
corresponding to the orbital period. The constant offset 
of the disturbance torques is attributable to the asym- 
15 metric orientation of solar panel 14 of FIG. la. In FIG. 
2, zero seconds and each 6,000 second increment there- 
after represent the times at which the spacecraft crosses 
the Earth’s Equator, known as the ascending node. 

FIGS. 3a, 36 and 3c represent as time functions the 
20 roll, pitch and yaw wheel momentum, respectively, in 
the absence of unloading. The changes in wheel mo- 
mentum with time are attributable to the disturbance 
torques illustrated in FIG. 2. Some components of the 
disturbance torques of FIG. 2 contribute to the cyclical 
25 variations of wheel momentum in FIG. 3, while others 
result in secular or continually increasing or decreasing 
components. For example, in FIG. 3a, the peak values 
of the cyclical component of wheel momentum in- 
creases with increasing time. At some point, therefore, 
30 such a continuing increase in amplitude will result in 
reaching the wheel momentum limit. Similarly, the 
pitch wheel momentum of FIG. 36 increases monotoni- 
cally with increasing time. In order to maintain continu- 
ous attitude control, the secular component of the 
35 wheel momentum must be unloaded. 

FIGS. 4a, 46 and 4c represent the roll, pitch and yaw 
torquer magnetic moment measured in ampere-tums- 
meter 2 (ATM 2 ) for a spacecraft using a prior art torquer 
control system which pulse width or continuously mod- 
40 ulates the torquer current. The prior art control system 
measures the Earth’s geomagnetic field B and the net 
momentum of the reaction wheel, and forms a vector 
cross product HXB. The normalized HXB vector is 
multiplied by the maximum magnetic moment which 
45 can be achieved along each of the roll, pitch and yaw 
axes to determine the torquer dipole command. In 
FIGS. 4a, 46 and 4c, large excursions or jitter in the roll 
and pitch moments, and lesser excursions in the yaw 
moments, occur at about 7500 seconds, 14,000 seconds, 
50 and 20,000 seconds, corresponding to passage of the 
spacecraft across the Earth’s magnetic poles. 

FIG. 5 is a simplified block diagram of a control 
system according to the invention. In FIG. 5, reaction 
wheel electronics 510 measure reaction wheel speeds ft. 
55 The ft signals are coupled to a block 512. A solar array 
drive is illustrated as a block 516, and the angle 0 of the 
solar array is coupled to a further block 518 which 
represents calculation of a momentum command vector 
h c . The momentum command vector is applied to block 
60 512, together with wheel speeds ft. Block 512 calculates 
momentum error h e . The momentum error h* calculated 
in block 512 is applied to an input of a block 514. A set 
of magnetometers for measuring three-axis geomagnetic 
field B is illustrated as a block 542. Block 542 is coupled 
65 to a further input of block 514. 

Block 514 represents jitter suppression and calcula- 
tion of a cross-product to produce a gain-controlled 
torquer magnetic moment command. The command is 
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applied to a further block 520, which represents effi- 
ciency zone and hysteresis logic. 

FIG. 6 is a simplified flow chart illustrating the calcu- 
lation of momentum command vector h c in block 518 of 
FIG. 5. In FIG. 6, the logic flow begins at a START 5 
block 610 and continues to a block 612, which repre- 
sents reading of the orbit rate vector co. The orbit rate is 
known from orbit dimensions, and may be prepro- 
grammed into a satellite or provided by an uplink from 
an Earth station. Block 614 represents the reading of the io 
solar array angle 0. Block 616 represents calculation of 
the spacecraft body momentum h&, which is the mo- 
mentum of the spacecraft body, excluding the solar 
array. The calculation is performed as a product, where 
I represents the inertia matrix of the spacecraft body 15 

h b =Iw (l) 
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is a vector representing the local geomagnetic fields. 
The logic proceeds to a block 814, which represents the 
calculation of the magnitude of vector B by taking the 
square root of the sum of the squares of the components 
along the X, Y and Z axes 


\B\ = W + By 2 + B; 2 


Block 816 represents the reading of momentum error h e 
from FIG. 7. Block 818 represents the calculation of 
| h e | by taking the square root of the sum of the squares 
of the components along the X, Y and Z axes 



Block 618 represents the calculation of the momentum 
ha of the solar array according to the equation 20 

h a =rXm((a>Xr) + r) (2) 

where 

r is the position vector from the spacecraft center of ^ 
mass to the solar array center of mass; 
x is the vector cross product; 
m is the mass of the solar array; and 
r is the time derivative of r, and both r and r may be 
functions of 0. ^ 

Block 620 represents the calculation of momentum 
command vector h c by summing 

hc = ho-\-hb (3) 

Finally, block 622 represents the storage of momentum 35 
command vector h c in order to make it available for 
further calculation. 

FIG. 7 represents the calculation of momentum error 
h e in block 512 of FIG. 5. In FIG. 7, the logic begins 
with a START block 710, and proceeds to a block 712 40 
which reads the value of fl for each of the N wheels of 
the reaction wheel assembly. The logic proceeds to a 
block 714, which represents the calculation of momen- 
tum for each of the wheels by multiplying the speed o> n 
of each wheel by the inertia I n for the corresponding 45 
wheel, to produce h 

— ’ (4) 

Block 716 represents the fetching or reading of momen- 50 
turn command vector h c from block 622 of FIG. 6. 
Block 718 represents the transformation of the momenta 
of the n reaction wheels to the spacecraft body frame. 
The transformation matrix is based upon the geometry 
of the various wheel axes, and may be stored as data in 55 
the spacecraft or transmitted to the spacecraft from an 
earth station. The transformed wheel momenta hwr otal 
is added to the momentum command vector he to pro- 
duce the desired momentum error h e in block 720: 

60 

ht=hc + h*TorrAL (5) 

Block 722 represents the storage of the calculated Value 
of h e to make it available for further calculations. 

FIG. 8 is a simplified block diagram representing the 65 
calculations performed in block 514 of FIG. 5. In FIG. 

8, the logic begins at a START block 810, and proceeds 
to a block 812, which represents the reading of B, which 


The logic then arrives at a decision block 820, in Which 
the magnitude of B from block 814 is compared with a 
limit value of B// m . If the absolute magnitude of the 
geomagnetic fields |B| is less than the limit value B/j m , 
the logic exits block 820 by the NO path and proceeds 
to a plurality of logic blocks illustrated together as 824. 
If | B | is greater than B// m , the logic flows from decision 
block 820 by the YES block to a further decision block 
822. Block 822 represents the comparison of | h e | with a 
limiting value h// m . If | h* | is greater than h// m , the logic 
leaves decision block 822 by the NO path and flows to 
the logic represented by block 824. Thus, so long as the 
absolute value of the geomagnetic field is less than the 
limiting value, which might be the case far from the 
Earth’s magnetic poles, or if the absolute value of mo- 
mentum error h e is greater than a limiting value, the 
logic flows by path 823 to the logic within block 824. 
Within block 824 the logic flows to a further block 828 
which represents calculation of an intermediate momen- 
tum value Mi as the quotient of vector cross products 

M\=heXB/\heXB\ (8) 

Logic block 830 calculates a further momentum prod- 
uct designated M2 

M 2 =m){M mQX ) (9) 

and this value of momentum is made available as M3 on 
a logic path 850 for further calculations. 

As mentioned, the magnitude of the geomagnetic or 
B fields near the poles is liable to be large while the 
magnitude of the net momentum is small. Under this 
condition, the logic flow through decision blocks 820 
and 822 of FIG. 8 will result in a logic flow from the 
YES output of decision block 822 to a block 836. Block 
836 represents the calculation of a system gain accord- 
ing to the Equation 

Gain=\h e \Bi im /h lim \B\ (10) 

The logic flows from block 836 to further logic con- 
tained within dotted lines 834. Logic blocks 838 and 840 
contained within logic block 834 are identical to logic 
blocks 828 and 830 contained within logic block 824. 
Product M2 produced at the output of logic block 840 is 
applied to a further logic block 842, which multiplies 
M2 by the gain calculated in block 836. The output of 
logic block 842 is coupled as M3 onto logic path 850. 
The values of M3 appearing on logic path 850 are alter- 
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natives which depend upon conditions. If the magnitude 
of the environmental magnetic field is small, or if the 
momentum error is large, M 3 has a value established by 
the processing of block 824, whereas if the momentum 
error is small and the magnetic field is large, the value 5 
is additionally multiplied by a gain function which is 
calculated in block 836. 

As so far described, the control system measures the 
magnitudes of the environmental magnetic field and the 
magnitude of the momentum error h e for controlling 10 
system gain. This prevents the jitter illustrated in FIG. 

4 , and provides improved pointing accuracy. 

FIG. 9 is a simplified flow chart of processing ac- 
cording to the invention. The arrangement of FIG. 9 a 
implements a hysteresis control which eliminates rapid j 5 
cycling at the transition between torquing and non-tor- 
quin g regions. Referring to FIG. 96, a momentum di- 
rection H' makes an angle of Blim\ with the measured 
environmental magnetic field B. The magnitude of 
angle Bum\ is selected so that the magnetic torquing is 2 q 
sufficiently effective or efficient. A suitable angle might 
be 45 degrees. If torquing were simply proscribed for 
angles less than Bum\ and permitted for angles greater 
than or equal to Olimu noise and/or measurement er- 
rors might cause repeated ON-OFF cycling of the 25 
torquers when the angle between h e and B is near 
Blimv This cycling is avoided by including a “dead 
zone” in the control. In FIG. 9b, the “dead zone” is the 
region between Blimx and Bum- The torquers are not 
actuated when the angle lies between Blimx and Bum . 30 

The logic of FIG. 9 a performs the operation de- 
scribed in conjunction with FIG. 9b, but for simplicity 
and computational convenience, the processing is per- 
formed in terms of parameters cf> = cos 0 , <J>£./a/i=cos 0 - 
Umx, and $um = cosBum - In FIG. 9b, the logic be- 35 
gins at a START block 910 and flows to a block 912, 
which represents setting of a flag T flag to a value of 
zero. T flag may assume logic values of zero or one, 
with logic zero indicating that the magnetic torquers 
were not ON during the previous iteration, and logic 
one indicating that they were previously ON. The logic 
flows through blocks 910 and 912 only during initial 
start-up. The logic flows from block 912 to a block 914, 
representing the calculation of as a quotient of prod- 


<j>=Bh e /\h e \-\B\ (ll) 

From block 914, the logic flows to a decision block 
916, in which the condition of T flag is examined. If 
TFL^Ghas a value of one (the torquers were previously 50 
ON), the logic exits from decision block 916 by the YES 
output, and flows to a further decision block 918. Deci- 
sion block 918 compares the absolute value of <J> with 
<$>UM\ If |<f> | is greater than §lim\, the logic exits from 
decision block 918 by the YES output. This corre- 55 
sponds to a condition in which B^Blimx in FIG. 9b. 
From the YES output of decision block 918, the logic 
flows to block 920, representing the setting of M com- 
mand to zero, which corresponds to a no-torquing com- 
mand. The state of Tfl^gIs set to zero in block 922, and 60 
the logic returns by a path 924 to block 810 of FIG. 8 in 
preparation for another iteration. 

If the torquers were previously ON, as described 
above, and <f> is less than, or equal to <t>i/Mi (B > Bum in 
FIG. 9b), the logic leaves decision block 918 of FIG. 9 a 65 
by the NO output, and flows to a block 934, represent- 
ing the setting of Mcommand to equal M 3 . This causes 
the torquers to operate at a magnitude established by 


8 

M 3 , where the magnitude may be accomplished by an 
appropriate percentage of pulse width modulation. 
T flag is set to a logic one in block 936, and the logic 
returns by path 924 to block 810 for another iteration. 

If the torquers were previously OFF the logic flow of 
FIG. 9a leaves decision block 916 by the NO output, 
and flows by a path 926 to a decision block 928. Deci- 
sion block 928 compares the absolute value of <f> with 
4>um ■ If 1 4> I < <t>Lim (corresponding to | B \ > Bum in 
FIG. 96) the logic leaves decision block 928 by the YES 
output, and flows by a logic path 930 to block 934, 
representing the setting of Mcommand to equal the 
value of M 3 . This results in torquing of the spacecraft. 
The logic then flows to block 936, representing setting 
of T flag to the ONE state. The logic returns by logic 
path 924 to block 810 to begin another iteration. 

If the torquers were previously OFF as mentioned 
above, the logic reaches decision block 928 of FIG. 9a. 
If | <j> | =§um (corresponding to B ^ Bum in FIG. 96), 
the logic leaves decision block 928 by the NO output, 
and proceeds by a logic path 938 to pass in sequence 
through blocks 920 and 922, representing no torquing 
and the setting Tflag— 0, respectively, following 
which another iteration begins. 

The dead zone operates in the following manner. If 0 
lies between Blim\ and Bum , B is greater than 0 LIM 1 
and less than Bum, corresponding to §<§lim\ and 
4> > 4 >l/M 2* The logic must flow through either decision 
block 918 or 928 in order to reach the torquing control 
condition, represented by blocks 934 and 936. When the 
logic flows through decision block 916 and §<<$>lim\, 
the logic always exits by the NO path, thereby flowing 
to blocks 934 and 936. When the logic flows through 
decision block 928, with (f> > 4>um, it always exits by 
the NO path to blocks 920 and 922. Since one path 
represents torquing and the other path represents no 
torquing, a dead zone exists. 

Other embodiments of the invention will be apparent 
to those skilled in the art. For example, the reaction or 
momentum wheels may exceed three in number, and 
may be in arbitrary orientations. The magnetic torquers 
may have any arbitrary orientation relative to the space- 
craft body, and may exceed three in number. Also, the 
magnetic torquers may have different sizes and/or dif- 
ferent magnetic moments. While the satellite orbit has 
been described as a polar orbit, the advantages of the 
system may be realized for other orbits. 

In addition, the torquing logic of FIG. 9 a can be 
implemented by weighting the commanded dipole, M 3 , 
by a factor, T, calculated as T = 1 — | <J> | where <j> is 
calculated as in block 914. 

What is claimed is: 

1. A method for maintaining satellite attitude, com- 
prising the steps of: 

measuring the direction of the ambient magnetic field 
in the vicinity of said satellite; 

determining the direction of the net wheel momen- 
tum of said satellite; 

determining the angle between said direction of said 
magnetic field and said direction of said momen- 
tum; 

comparing said angle with a predetermined value; 
and 

energizing magnetic torquers for torquing said satel- 
lite for reducing said wheel momentum only when 
said angle exceeds said predetermined value. 
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2. A method according to claim 1 further comprising 
the step of measuring the magnitude of said ambient 
magnetic field; and 

in said step of energizing magnetic torquers, perform- 
ing the further step of controlling the magnitude of 5 
the energization in response to said magnitude of 
said ambient magnetic field. 

3. A method for maintaining satellite attitude, com- 
prising the steps of: 

measuring the direction of the ambient magnetic field 10 
in the vicinity of said satellite; 

determining the absolute magnitude of said ambient 
magnetic field; 

determining the direction of the net wheel momen- 
tum of said satellite; 15 

determining the angle between said direction of said 
magnetic field and said direction of said net wheel 
momentum; 

comparing said angle with a predetermined value; 

establishing the magnitude by which magnetic torqu- 20 
ers are to be energized in response to said magni- 
tude of said ambient magnetic field, by the further 
steps of: (a) determining the absolute magnitude of 
the sum of the spacecraft body, ancillary panel and 
wheel momenta to generate a momentum error 25 
magnitude; (b) comparing said absolute magnitude 
of said ambient magnetic field with a limiting value 
of said field; (c) comparing said momentum error 
magnitude with a limiting momentum magnitude; 


10 

and (d) multiplying said magnitude of said energi- 
zation by a multiplier other than unity only when 
said absolute magnitude of said ambient magnetic 
field exceeds said limiting value of said field and 
said momentum error magnitude is less than said 
limiting momentum value; and 

energizing said magnetic torquers for torquing said 
satellite for reducing said wheel momentum only 
when said angle exceeds said predetermined value. 

4 . A method according to claim 3 further comprising 
the step of: 

at least when said absolute value of said ambient 
magnetic field is less than said limiting value of said 
field and said momentum error magnitude is 
greater than said limiting momentum value, multi- 
plying said magnitude of said energization by a 
value of unity. 

5. A method for maintaining satellite attitude, com- 
prising the steps of: 

measuring the direction of the ambient magnetic field 
in the vicinity of said satellite; 

determining the direction of the net wheel momen- 
tum of said satellite; 

determining the angle between said direct of said 
magnetic field and said direction of said momen- 
tum; and 

torquing to reduce wheel momentum in proportion to 
the magnitude of said angle, 
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